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the origin, etc... which were obtained since 1977 in the framework of the Schrodinger equation 
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1 Introduction 



Everybody knows what was called, in particle physics, the "October revolution" because it 
took place in November 1974: the discovery of the J/ip particle PQ, which was rather quickly 
unerstood to be a quark-antiquark system, cc, where c, the charmed quark, was precisely the 
quark predicted in 1970 by Glashow, Iliopoulos and Maiani [3]. It was proposed to treat the 
system as a non-relativistic system satisfying the Schrodinger equation because the charm quark 
was heavy, 1 to 2 GeV/c 2 . Whether this was really allowed or not will not be discussed here. 
The fact is that, as we shall see later, it was very successful. 

In January 1976, I was invited to attend the "Orbis Scientiae" conference in Coral Gables 
where I was planning to give a talk on the diffractive peak at high energy, and on my way, I 
stopped for a few days at Rockefeller University. There, the late Baqi Beg told me: "You are an 
expert on the Schrodinger equation. Can't you explain why all models of Charmonium predict 
that in between the ground state J/ip and what is believed to be the first radial excitation, the 
ip', there should exist a multiplet of P states". 

I thought about this problem and it is only in 1977 that I found an imperfect answer 
to this question, namely a reasonable condition on the quark-antiquark potential which would 
guarantee this above-mentioned property. Naturally the existing models satisfied this condition. 
It should also be said that the P states were discovered at the right place. It is only much later, 
in 1984, that with Baumgartner and Grosse, I found a perfectly natural and beautiful condition 
to ensure the correct order of energy levels. But let us come back to the Spring of 1977. I had 
published one paper on the level ordering [3] and a second one on the relative magnitudes of the 
wave functions at the origin of the J/ip and ip' [I]- Professor A. Zichichi had his office exactly 
opposite to mine. He read my papers, and thought they were interesting. Then he asked me to 
give several lectures at the School of Subnuclear Physics in Erice on all aspects of Charmonium 
during the summer of 1977. I was completely unprepared to do that, and my wife is a witness 
that I was completely panicked. The summer came, and the lectures went relatively well [5]. 
In the end it turned out that being forced to give these lectures was a blessing in disguise, a 
blessing of God, God being personified by Antonino Zichichi! Indeed in the following years, an 
important fraction of my activity was devoted to Charmonium, or rather Quarkonium, since bb 
systems had also been discovered, and since ss, cs, bs, could also be meaningly included among 
heavy quark-antiquark systems. Not only mesons but also baryons containing heavy quarks 
were studied. 

It is very difficult, both because of the abundance of material and because of the complexity 
of certain proofs to present all the results here. I have already said that the theorems proved 
in 1977 were superseded by the much nicer results of 1984, and so, we will forget about these 
early theorems. 

I would like to explain that later my activity or rather "our" activity, because I had many 
collaborators, was divided into pursuing the derivation of rigorous results on the energy levels 
and wave functions of systems satisfying the Schrodinger equation (with possible applications 
outside quarkonium physics, like atomic physics, muonic atoms), and of phenomenological fits 
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of the quarkonium spectra which have happened to possess a very impressive predictive power. 



Among my collaborators I would list, in chronological order of appearance, Harald Grosse, 
Reinhold Bertlmann, Jean-Marc Richard, Alan Common, Bernard Baumgartner, Joachim 
Stubbe, Jon Rosner. 

Here I want to summarize both aspects of our activity. Concerning rigourous results, there 
has been already a Physics Reports published by Grosse and myself [5] following a Physics 
Reports by Quigg and Rosner [7] on the same subject. There was also a review presented 
at the School of Subnuclear Physics in Erice in 1992 [Sj. Later, Grosse and myself produced 
a fairly complete book j9] recently reprinted in paperback. This book contains results which 
concern not only quarkonium but also muonic atoms and alcaline atoms to which our theorems 
apply. 

In the last section we present an update of the fits of heavy quark systems (quarkonium and 
heavy baryons) with a simple model. This is an improved version, containing new experimental 
data, of a review at the Montpellier conference in 1997 [TO] . 

2 Order of Energy Levels 

Figures 1 and 2 present the energy levels of the cc and bb systems respectively. One sees 

i) that the £ = 1 states are between the £ = states, 

ii) that the average energy of the £ — 1 states is larger than the average of the energies of 
the £ = states immediately above and below, and that the £ = 2 state of the cc system 
is above the first radial excitation. 

What was found to be the "good" condition to explain i) is that the central potential should 
be such that its Laplacian is positive [IT] . 

We can state the following theorem: 

Theorem I 

Let E(n, £) be an energy level in a central potential, n being the number of nodes of the radial 
wave function, £ the orbital angular moment. Then 

E(n+1,£)> E(n,£ + 1) if r 2 A V(r)> . (1) 

Now: why is this a good criteria, for the quark-antiquark potential. The answer is that this is 
a strong version of asymptotic freedom. 

Call the force between a quark and an antiquark 
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In the Coulomb case, a(r) is essentially constant and 



i.e., AV^ = outside the origin. 

In QCD, asymptotic freedom can be expressed as 

j- <*{r) > , (4) 
dr 

where a(r) is the running coupling constant and hence AV(r) > 0. 

Now let us come to ii). We know, since Newton (according to Markus Fierz), that there 
are two and only two kinds of forces for which all bound states have classical closed orbits: the 
Coulomb and harmonic oscillator forces. 

Theorem I represents a comparison of the potential with the Coulomb potential. It might 
be interesting to compare the quark- ant iquark potential with the harmonic oscillator potential, 
V = r 2 . V = a + br 2 satisfies 

d 1 dV _ 
dr r dr 

One can prove the following theorem [12]. 

Theorem II 
If 

d_ 1 dV >Q 

dv v dv 

E(n+1,£) >E(nJ + 2) (5) 

In the case of quarkonium, it has been shown by E. Seiler [13], from lattice QCD, that the 
quark-antiquark potential is concave and increasing. Hence 

d 1 dV 1 dV 1 d 2 V 
Tr r ^ = ~~ 2 -d7 + r ^ <0 

In this way, we understand that the mass of the if) (n = 0, i = 2, ) is higher than the mass 
of the ip' (n = 1, i = 0). We shall come back later on the question of why the I = 1 states are 
above the average of the neighbouring two £ = states. 

Outside these two theorems on the level ordering there are others, where the criterium is 
the sign of D a V(r) where 

D a = ^ + (5-3a)- -j- + 2(1 — a) (2 - a) \ (6) 

ar r dr r z 

the previous cases correspond to 

a = 1 and a = 2 
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These theorems are not optimafj], but we have conjectured with a relatively strong basis that if 



r-£-(« 2 -3)-f| V(r)<0 (7) 



dr 2 dr 
with a > 2, then £(n, < £(n — 1, £ + a), for instance if we take 

V(r) = r 4 

£(n,£) < £(n- + ^6) 

Indeed, for instance 

£(1,0) =11.6, £(0,\/6) =12.6 
£(1,5) =35,7, £(0,5 + ^) =36.3 

In another special case, V(r) = r 5 a rigorous but painful proof gives 

E(n = 0,£ + V7) > E(n = 1,1) . 



Let us return now to Theorem I. So far it was applied to quarkonium, but there are other 
applications. The first one is muonic atoms. In muonic atoms, the Bohr radius is so small that 
the muon "sees" the extension of the nucleus. 

Since the nucleus, in the first approximation where it is supposed to be made of protons and 
neutrons, has a positive charge distribution, the potential exerted on the muon has a positive 
Laplacian and 

E(n + 1,£) > E(n,£+1) 

This is indeed what was found in early calculations of muonic energy levels [H]. However, 
this was in fact independent of the details of the chosen charge distribution. This is seen 
experimentally. For instance, using the standard atomic spectroscopic notation, N, £, total 
angular momentum, with iV = n + £ + 1, for 138 Ba) atoms [15] 

2s 1/2 - 2p 1/2 = +405.41 KeV 
3p 3 /2-3d 3 /2 = 8.69 KeV 



One might object that relativistic corrections cannot be neglected. However, it has been 
shown [16J (at least for perturbations around the Coulomb potential) that the inequalities 
survive for the Dirac equation as long as one compares levels with the same total angular 
momentum. 

The other field where Theorem I is useful is atomic physics, in the approximation where 
the alcaline atoms can be treated as one outer electron plus a closed shell. The nucleus, seen 
from far away by the outer electron is now essentially pointlike and the effective potential is 
a Coulomb potential plus the result of a negative charge distribution. So it has a negative 

1 For details, see [5], p. 43 
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Laplacian. Take for instance Lithium. It has three electrons. Two of them form a closed shell. 
What will the outer electron be: 

2s or 2p ? 

It will be 2s because of the negativity of the Laplacian. In this way, Theorem I helps to 
understand many aspects of the Mendeleiff classification which, previously, were explained 
by not very convincing handwaving arguments. Figure 3, taken from [T7j shows that these 
considerations hold also for excited states of the lithium atom. 
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One sees also that for large principal quantum numbers the Coulomb degeneracy tends to 
reappear. On Fig. 3, it is interesting to see that 

3s < 3p < 3d 

more precisely, one has, theoretically [18], and experimentally [T9] : 

3si /2 < 3pi /a < 3p 3/2 < 3d 3/2 < 3d 5/2 

27206.066 < 30925.517 < 30925.613 < 31283.018 < 31283.053 cnT 1 (8) 

These inequalities are valid for an outer electron satisfying the Dirac equation with a 
monotonous increasing potential such that AV < 0. The assumption of monotonicity is prob- 
ably superfluous. 



3 Spacing of Energy Levels 

Since energy levels depend on two quantum numbers, n and £, the problem is very broad and 
very difficult to explore completely. We shall restrict ourselves to two extreme cases: 

a) Compare the spacing of the energy levels for fixed n, and in fact we shall restrict ourselves 
to n = 0, i.e., the purely angular excitations. 

b) Study the relative spacings of the energy levels for fixed £, specifically £ = 0. 

While, for case a) there is a large amount of exact results, for case b) we have only imperfect 
indications. The results obtained in a), however, will allow us to answer question ii) in Section II. 

In these problems a simple reference potential is the harmonic oscillator V = r 2 . There all 
the fixed n, increasing £ spacings are equal and all the fixed £ increasing n spacings are equal. 

a) Spacing of purely angular excitations n = 0, increasing £. There one gets easily 

Theorem HI %2U\ 
If 

d 1 dV 
dr r dr 
E(0,l+2)-E(0,l+l) 
E(0,£ + 1)-E(0,e) {) 

There exists also a differential form of this inequality. 

If we combine this with Theorem II, which gives, under the same conditions 

E(l,£) < E(0,£ + 2) 
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we get 

E{0,£+1)>~ [E(0,£) + E(1,£)} 

and the special case 

E(0,1)>~[E(0,0) + E(1,0)} 

which is precisely what is observed in the cc and bb spectra. 

We would also like to have the analogue inequality for the second £ = 1 state, E(l,l) but 
so far this is not the case. 

A more subtle inequality for angular excitations is obtained if AV(r)< 0. Denote E(0,£) as 
Ee, then 

Theorem IV [21] 
If 

d 9 dV ^ 
— r 

dr dr < 



E i+l - E e > 21 + 3 / £ ' 
E t -Et-i < 2£+ 1 \£ + 2 



(10) 



where the right-hand-side corresponds to the pure Coulomb case. Application to the cc system 
is rather frustrating. We get 

E 2 -E 1 ^ 5_ 
E\ — Eq 27 

trivially satisfied by experiment. 

Muonic atoms, on the other hand, are much more interesting. They also have r 2 AV > 
and the spacing between purely angular excitations is what is most easily observed by looking 
at cascading 7 rays by favoured electric dipole transitions. For small nuclear charge Z or for 
large angular momenta the deviations for the Coulomb value are small while they are large for 
small I and large Z. 



After relativistic corrections, we getd 



Tabic 1 



z 




E4f — Esd 

Esd-E2 v 


E 4 f-E 3 d 


z - 


♦ 0.185 


0.350 


0.463 


20 


0.202 


0.350 


0.463 


40 


0.253 


0.350 


0.463 


60 


0.333 


0.353 


0.463 


80 


0.462 


0.370 


0.463 



2 The references on the experimental material used in Tables 1, 2, 3 and 4 can be found in Ref. [21] . 
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Conversely, for Alcaline atoms we have AV < 0. Here we just give the Lithium sequence, 
i.e., ions with three electrons and increasing charge. Here, it is just the reverse. When Z 
gets large, the Coulomb limit is approched because the charge of the nucleus dominates on the 
electron cloud. For extremely large Z (this is fictitious) the orbits are inside the cloud and do 
not see it. 



Tabic 2 



z 


E3d — E2p 




3 


0.326 


0.462 


6 


0.329 


0.462 


9 


0.334 


0.462 


12 


0.338 


0.462 


15 


0.339 




oo 


0.350 


0.463 



Finally we give inequalities on the spacing between states with the same quantum numbers 
and different total angular momentum, i.e., the fine splitting. If, keeping n = 

8 = E(£,J = £+±)-E(l,f = l-±) (11) 



s(e + iy i(i + i) 3 
W) < (^ + 2) 4 ' " 



(12) 



we have in the semi-relativistic approximation 

§ _ 21 +1 ldV\ 
Am? \dr I 

where the expectation value is to be taken using the Schrodinger wave function. 

We have two theorems [21~] . 
Theorem V 

'M^W = Si {E ^- E ' f (13) 



and 

a// ii /// ■ rvi 

if r 2 AV(r)> (14) 



< 

for muonic atoms , i.e., r 2 AV(r) > 0, we give only a sample of the results and forget the 
uncertainties to improve legibility 



Tabic 3 



z 


5(1) 


A(l) 


8(2) 


fi > °- 098 




in keV 






26 


4.2 


4.4 


0.47 


0.112 


33 


11.10 


11.45 


2.00 


0.18 


41 


23.15 


29.10 


2.68 


0.116 


50 


45.7 


60.60 


5.65 


0.123 
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Applications to Alcaline atoms, with AV < 0, are on the other hand rather disappointing. 
Let us remember that Eq. (12) is violated by the levels of the Sodium atom since the sign 
of the spacing of the well-known Sodium doublet (which produces by transition to the ground 
state the horrible yellow light) is the opposite. Exact treatment by the Dirac equation does not 
help. It is a typically many-body effect. 

Concerning the lithium isoelectronic sequence, i.e., ions with three electrons, the situation 
is not as bad in the sense that 5(1) and 5(2) have the right sign. However, the inequality 

5(1) > A(l) and 5(2 > A(2) (15) 

are violated for Lithium itself. It is only for Z > 6 (Be II, etc.) that they are satisfied. This is 
seen in Table 4 (the units are cm -1 ) 



Table 4 



z 


5(1) 


A(l) 


5(2) 


A(2) 


d M < °- 098 




in cm -1 








3 


0.337 


0.422 


0.037 


0.036 


0.1097 


5 


34.1 


34.2 


3.1 


2.5 


0.0909 


6 


107.1 


106.6 


10.5 


9.278 


0.098 


9 


975.8 


962.7 


90.0 


87 


0.092 


12 


3978 


3938 


470 


362 


0.0118 


15 


11310 


11100 


1000 




0.088 


24 


90910 


90000 









We see also that 5(1)/ A(l) is always very close to unity. There are also inequalities of the 
same kind for hyperfine splittings but, for these, we refer to the original work [21J. 

We turn now to the problem of spacings between £ = levels [22] • We present a number of 
incomplete results, which give strong indications, but completely rigorous proofs are lacking. 
Two things are certainly true: 

a) the £ = levels of the harmonic oscillator potential are equally spaced, as everybody 
knows; 

b) the £ = levels of the linear potential have a spacing decreasing with n. 

We study now the neighbourhood of a), i.e., potentials close to the harmonic oscillator 
potential. 

A preliminary remark is that the spacing between the n = levels will remain constant 
not only if we replace r 2 by A r 2 + B, but also if we add a C/r 2 term. Adding such a term 
is equivalent to changing the angular momentum by an (generally not integer) amount £ such 
that C = £(£ + 1). Since the "Regge trajectories" of the harmonic oscillator are linear and 
parallel, the spacing between the energy levels will not change. 
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An interesting quantity to control the spacings is 

Z(V,r) — ^ r 5 ^ 1 ^ 

' <ir dr r dr 

Z vanishes for 

V = const., V = r 2 , V= \ . (17) 

What we have proved, using a new kind of raising and lowering operators is that if 

V = r 2 + \v , (18) 

for A sufficiently small , the spacing between the energy levels increases with n if Z(V, r) > 
Vr and the spacing between the energy levels decreases with n if Z(V, r) < Vr. 

However, in the non-perturbative case, we do not know what happens. In fact we even have 
counter examples. 

Take the "partially soluble" potential [23J 

V = r 6 - 9r 2 

the first five energy levels are given by 

E , E 4 = =f^/480 + 96vTT 

Ei, E 3 = =F-y/ 480 - 96^TT ( 19 ) 
E 2 =0 



Obviously, even though Z(V,r) is positive, the levels show a symmetric pattern around 
E2 = 0. However, there is no known counterexample in which V is monotonous increasing. 

We are tempted to make the conjecture that if Z > and if the potential is monotonous, 
the spacings increase, and if Z < and the potential is monotonous, the spacings decrease. 
"Experimental" tests with V = r a indicate that the spacing increases with n if a > 2 and 
decreases with n if a < 2. 

We study now the neighbourhood of b), the purely linear potential. Then the solutions of 
the Schrodinger equation are Airy functions. Then, the energy levels are given by E n = — r n _i, 
where r n is the n th zero of the Airy function. It is obvious that the spacing between r n and 
r n+ i is smaller than the spacing between r„_i and r n because the potential is stronger. 

We strongly believe that for all concave potentials, which correspond to the physical sit- 
uation for quark-antiquark systems [13] . the spacing between £ = levels decreases with n. 
This is precisely what is indicated by the WKB approximation, in which we can regard as a 
continuous variable. The spacing will decrease if 

d 2 E d 2 n 

< , i.e. — - > . 

d n 2 dE 2 
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Indeed, in the WKB approximation 

d 2 n 



dE 2 27T v(o)y/E-v(p) 

Ie-v>o v' 2 ve-v d r (20) 

which is positive if V'(0) > and V" < 0. 

The bb spectrum possesses this property, at least below the BB threshold |24j : 

My - M 7 = 560 MeV 
My - My = 332 MeV 

4 The Wave Function at the Origin, the Kinetic Energy, 
etc. 

Perhaps it is worth mentioning results on the wave function at the origin because one of the two 
papers which initially attracted the attention of Professor Zichichi is precisely on this subject 

m- 

The wave function at the origin appears in the so-called Van Royen-Weisskopf formula [25] 
(also proposed by Pietschmann and Thirring, M. Krammer and H. Krasemann) which controls 
the leptonic width of quarkonium: 

T e+e - = IGvr oc e% ME , (21) 

where cq designate the charge of the quarks in units of e. |f/>(0)| 2 also controls the hadronic 
width, i.e., the decay into three gluons of quarkonium. The reduced wave function at the origin 
is given for the I = wave function by a formula attributed to Schwinger: 

o r°° dV 

K(0)) 2 = 2mj o — ul dr , (22) 

where 

rib 

u - 



,'Att 

with 

ul dr = 1 



if V is linear , = const, and hence, because of the normalization of u n , \u' n (0)\ and |^ n (0)| 
are independent of n. 

We proved the following theorem: if V is concave, i.e., if V < 

hM0)| 2 <hM0)| 2 , (23) 
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if V is convex, we get the reverse. 



Indeed, even if you take into account the change of mass in the Weisskopf-Van Royen 
formula, experiment indicates that the leptonic widths of the J/ip and ip' fit with a concave 
potential, which is what we expect from lattice QCD [13J. 

Can we go beyond that result? All we can say is that we have strong indications that the 
leptonic widths decrease with n if V is concave. For instance we can prove that \ip n {0) | 2 goes 
to zero for n — > oo if V is concave. 

Outside the wave function at the origin there are other quantities of interest on which we 
have a control, such as the mean kinetic energy, the root mean square radius, the electric dipole 
transition matrix elements, etc. For all these we send the reader to Ref. [9]. Let us just mention 
that this allows to set constraints on the "Schrodinger mass" of the quarks, if the potential 
is flavour independent, from the values of the I = and I = 1 energy levels of the cc and bb 
systems: 

m b -m c > 3.29 GeV , (24) 

while a naive approach gives 

m b - m c > - {M bl - M cZ ) = 3.18 GeV (25) 

The mean kinetic energy of the system, T, that we call T(m), where m is the common quark 
mass in quarkonium, appears in the change in binding energy of a quark- ant iquark system when 
the mass changes for a fixed potential (called flavour- independent for quarkonium), which is, 
from the Feynman-Hellman theorem 

d . T(m) , , 

— E(m) = — 26 

dm m 

This proves already (25), from the positivity of T. 
Now, with M > m 

rM T(JX) 

H 

From the concavity with respect to any parameter entering linearly in the Hamiltonian, we get 
that mT(m) is increasing with m. In fact, for a square well potential mT(m) is constant. 



E(M) = E{m) - f M dfi 



However, if the potential is concave, which we believe from lattice QCD [13], we get a 
stronger result [26] : 

m 1/3 T(m)/ (27) 

with m. 

On the other hand, an inequality on T(m) has been obtained [27] 

T(m) > - [E(n = 0,£=l,m)-E{n = 0,£ = 0,m] , (28) 
which is saturatedby a harmonic oscillator potential. 
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From these considerations it is possible to get a lower limit to the mass difference between 
the b quark mass and the c quark mass, taking into account that the mass of a QQ system is 
2m,Q + E{mQ): 

m b -m c > 3.32 MeV (29) 

Without (27), we would get the weaker result (24). This is only a sample of the many results 
obtained in this domain, for which we send the reader to Ref. [9]. 



5 A Fit of Heavy Quark Systems by a Potential Model 

Now, we are leaving mathematical physics and turning to (dirty) phenomenology. After the 
discovery of the Upsilon and Upsilon prime systems by Lederman in 1977 [28J, Quigg and Rosner 
made the remark that the spacing of the levels of the cc and bb systems is almost the same. If 
a potential model is acceptable, and if this potential is flavour independent, it is tempting to 
say that the potential could just be logarithmic, i.e., C£n r. Then it is almost obvious that the 
spacings of all energy levels will be independent of the mass, because re-scaling r to take into 
account the change of the mass will shift the potential by a constant. In 1981, I tried a small 
generalization of this by taking V — A + Br a , and adjusting A, B and a to the known levels. 
I added too a Fermi-like hyperfine splitting to "explain" the J/tp — r] c mass difference. This 
seemed to be successful [29] , [30] , for mesons made of b, c and even s quarks and antiquarks. 
At the same time, Jean-Marc Richard, using the "rule" Vqq = | Vqq reproduced beautifully 
the mass of the Q~ baryon |31j . 

Then, my model was taken with disbelief by many physicists, including some friends, es- 
pecially Russians, who could not understand how a naive potential model could be justified. 
Most people preferred QCD sum rules, or lattice QCD, still in infancy. 

However, I took the position that it was worth experimenting with this model even if it 
was difficult to justify. With time, many new particles and energy levels were discovered and 
happened to fit perfectly with this model. I remember a seminar by Lorenzo Foa, announcing 
the discovery of the B s meson by Aleph at LEP and saying that it was not necessary to give its 
mass because I had already predicted it. With, may be, one exception, this continues till now. 

In fact, a first version including only b and c quarks was proposed [29]. However, seeing 
the incredible success of the fit, Murray Gell-Mann suggested that one should go further in the 
"absurd" (since c quarks are not really non-relativistic) and include the strange quark. Let 
me give the numerical elements of the fit, made in 1981 |30j . with the existing experimental 
information. The potential is 

V = A + Br a 

the fit gives 

A = -8.064 
B = 6.870 , 

in units which are powers of GeV 

a = 0.1 
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m b = 5.174 
m c = 1.8 
m s = 0.518 

(notice that this agrees with the lower limit (29) on mj — m c ) . 

Futhermore, the large hyperfine splitting between J/ip and r) c has to be taken into account 
by a phenomenological spin-spin interaction 



C 



m!m 2 

C is adjusted to the J /i/j — r\ c splitting. 



{a x .a 2 ) 5 (n - f 2 ) 



Overall this is a seven parameter fit. We do not try to predict the fine splitting between 
the i = 1 states and just predict the weighted average. 

Table 5 gives the results. It contains 30 experimental numbers. It is an updated version of a 
table presented at the Montpellier 1997 conference [10]. The figures with stars are experimental 
results which were not known in 1981 (ten stars!). The fact that the T]' c is 15 MeV higher than 
predicted has been understood long ago because of the vicinity of the DD threshold (by J.-M. 
Richard) . 



16 



Table 5 

Masses and relative leptonic widths for cc, bb, ss, 
cs, bs, be, and of the fl and fl c baryons. 



Quark 
System 


States 


Theory 


Experiment 


Theory 


Experiment 


cc 




3.095 


3.097 


1 


1 






3.687 


3.686 


0.40 


0.46 ± 0.06 




T 


4.032 


4.040 


0.25 


0.16 ± 0.02 




Xc average 


3.502 


3.525 








and h c 












ip" average 


3.787 


3.770 








Vc 


2.980 


2.980 








rfc 


3.641 


3.656* 






bb 


T 


9.46 


9.46 


1 


1 






10.02 


10.02 


0.35 


0.32 ± 0.05 




T" 


10.36 


10.35 


0.27 


0.24 ± 0.05 






10.60 


10.58 


0.21 


0.23 ± 0.06 






10.76 


10.86 * 








■y b average 

/V o 


9.86 


9.90 * 










10.24 


10.26 * 






ss 





1.02 


1.02 








fi average 

J J- o 


1.42 


1.44 








P state 












0' 


1.634 


1.680 






cs 




1.99 


1.97 










2.11 


2.11 * 








-D** average 


2.537 


2.536 * 






bs 


5, 


5.354 


5.369 ± 0.005 * 








B* s 


5.408 


5.416 ± 0.004 * 






be 


B c 


6.25 


6.4± 0.39 ± 0.13 * 










6.32 








n 




1.666 


1.672 










2.708 


2.697 ± 0.0025*) 






These 


experimental numbers have been obtained after the initial fit in 1981. 
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We turn now to the baryon sector. Here we have been using the rule Vqq = \ Vqq. Why 
this rule? A three quark system must be colour singlet. Hence in a three quark baryon, a quark 
pair must be a 3 colour state. So, if two quarks are close together, the third quark sees them 
as a 3 state, i.e., an antiquark. Hence the potential between the third quark and the quark 
pair is Vqq. Dividing by two we get a potential | Vqq between the third quark and each of the 
quarks of the pair. 

This may seem doubtful but worth trying. Jean Marc Richard has done this for the Q 
particle made of three strange quarks [31]. This gives Mq = 1666 MeV, to be compared with 
the experimental value 1672 MeV. Encouraged by this result we predict [32] 

M Qc = 2708 MeV, 

which is now, experimentally [33] 2697.5 ± 2.5 MeV, only 10 MeV away. 

Other predictions, still to be tested, are 

M Ucc = 3.737 MeV 
Mn ccc = 4.787 MeV 

The ccc system, according to Bjorken [34"] . is one of the most interesting quark systems. It 
is stable with respect to strong interactions and has a lifetime of the order of 2 x 10~ 13 seconds. 
It might be difficult to produce, but, who knows, it might be seen at the LHC. After all, at 
LEP, many interesting particles such as B s , B c or Af, were seen while their observation was not 
initially planned. 
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